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Abstract 

I derive analytically the spectrum of the CMB fluctuations. The final result for C\ 
is presented in terms of elementary functions with an explicit dependence on the basic 
cosmological parameters. This result is in a rather good agreement with CMBFAST for 
a wide range of parameters around concordance model. This allows us to understand the 
physical reasons for dependence of the particular features of the CMB spectrum on the 
basic cosmological parameters and to estimate the possible accuracy of their determination. 
I also analyse the degeneracy of the spectrum with respect to certain combinations of the 
cosmological parameters. 



1 Introduction 



After recombination the primordial radiation doesn't interact anymore with the matter and most 
of the photons come to us without further scattering. Since the radiation is extremely isotropic 
in nearly all angular scales we conclude that at the moment of recombination the universe was 
extremely homogeneous and its temperature could not vary from place to place more than about 
few times in thousandth of the percent. 

On the other hand, the origin of the large scale structure requires the presence of small 
inhomogeneities in the distribution of the matter and therefore the temperature of CMB should 
also vary a little bit. These variations are observed today as the angular fluctuations of the CMB 
temperature PQ. The expected fluctuations in a given angular scale are basically determined by 
the inhomogeneities on the spatial scales having today an appropriate angular size if placed at 
the distance corresponding to the recombination redshift. 
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The Hubble scale at recombination epoch plays especially important role, distinguishing 
the inhomogeneities which are still frozen from those ones which already entered the horizon and 
therefore could be amplified by gravitational instability. At the scales bigger than the Hubble 
size, the perturbations generated during inflation remain unchanged. Therefore, observing the 
fluctuations on the angular scales 6 > 1°, corresponding to super-Hubble scales at recombination, 
we directly probe the primordial inflationary spectrum not influenced by the evolution. The 
perturbations which entered the horizon before recombination evolve in a rather complicated 
way. The transfer functions relating the initial spectrum to the resulting one strongly depends on 
the major cosmological parameters, and the shape of the CMB fluctuations spectrum at 6 < 1° 
is very sensitive to the exact values of these parameters. Therefore by measuring the fluctuations 
at small angular scales we can determine these parameters. 

The recent observations of the CMB fluctuations P give us a hope that finally we will 
be able to determine the cosmological parameters with a very high precision. One of the most 
important parameters among them is the spectral index n s characterizing the initial perturba- 
tions. According to inflationary paradigm n s should deviate from n s = 1 and be in the range 
0.92 < n s < 0.97 depending on the particular scenario of the simple 1 inflation |2], [H|. It is very 
important to find these deviations to confirm or disprove inflationary paradigm. The accuracy of 
the current observations is not yet high enough to conclude about the deviations of the spectral 
index from the flat one jl]. However the future measurements will be able to reach the needed 
precision. 

The CMB spectrum depends on the various cosmological parameters in a rather com- 
plicated way. It is very important to clarify this dependence to be sure which features of the 
spectrum are most sensitive to the particular combinations of cosmological parameters. The 
usual approach using the computer code CMBFAST 5J is very helpful, but it does not com- 
pletely solve the problem since the parameter space has too many dimensions. There are various 
semi-analytical and analytical approaches to this problem |Sl,jZj- However I was not able to find 
in the literature elementary analytical expression which would explicitly describe the dependence 
of the CMB spectrum on the cosmological parameters and would be in a reasonably good agree- 
ment with numerics. In this paper I derive such expressions. The main results of this paper are 
the equations (|9*2>(ITUUJ) . 

I start with a pedagogical introduction reminding the derivation of the Sachs- Wolfe effect 
in the conformal Newtonian coordinate system |3] and first make the calculations assuming the 
instantaneous recombination. In this approximation the radiation can be well described in a per- 
fect fluid approximation before recombination and as an ensemble of free photons immediately 
after that. This is well justified by causality only when we consider the fluctuations corresponding 
to the superhorizon scales. At small angular scales the delayed recombination is quite important 
and leads to an extra damping of the fluctuations. Therefore as a next step I show how the for- 
mulae obtained in the instantaneous recombination approximation should be modified to account 
for this effect. Finally the spectrum for the small angular scales is derived and the precision of 

1 Under simple inflation I mean the scenarios with the minimal number of free parameters. 
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the determination of cosmological parameters and degeneracy of the spectrum with respect to 
certain combinations of these parameters is discussed. 



One important simplification I make is that I nearly always consider the most observa- 
tionally favored case of a flat universe. The modifications of the most important features of the 
CMB spectrum due to the spatial curvature are rather obvious. 

In Appendix A I derive the analytical formulae describing non-instantaneous recombina- 
tion which I use in the section on the finite thickness effects. In Appendix B the transfer functions 
in short and longwave limits are derived in the conformal Newtonian gauge. 



Before recombination the radiation is strongly coupled to the matter and it can be well described 
by a perfect fluid approximation. When the hydrogen becomes neutral, most of the photons 
do not interact anymore with the matter and therefore to describe them we need the kinetic 
equation. 

The free propagating photons are described by the distribution function / defined via 



where dN is the number of particles at time r] in the appropriate element of the phase volume 
d 3 xd 3 p = dx 1 dx 2 dx 3 dpidp2dp3, so that / is the particle density in the one-particle phase space. I 
assume that the indices a, (5... run always over 0, 3 while i, k take only spatial values 1, 2, 3. The 
phase volume is invariant under coordinate transformations and hence the distribution function / 
is a space-time scalar. Since the phase volume is conserved along the trajectory, the distribution 
function in the absence of the scatterings should obey the collisionless Boltzmann equation 



where dx l /drj and dpi/drj are the appropriate derivatives calculated on the photon's geodesic. 

Temperature and its transformation properties. The energy (frequency) of the photon 
with the 4-momentum p a measured by an observer having the 4-velocity u a is equal to the scalar 
product of these vectors: cu = p a u a . This can be easily understood by going to the local inertial 
coordinate frame of the observer. If the spectrum of the quanta coming to an observer from a 
particular direction on the sky, characterized by the vector ni = —pi/p, where p = (T.pf) 1 ^ 2 , is 
the Planckian one then the temperature, defined via 



2 Sachs- Wolfe effect 



dN = f(x\pj,r])d 3 xd 3 p 



(1) 




(2) 




(3) 
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generically depends not only on the direction n« but also on the moment of time r\ and the position 
of the observer x % . The factor two here accounts for two possible polarizations of the photons. 
From now on I consider the Universe where the fluctuations of the temperature are very small 
and therefore one can write 

T(x a ,k) = T ( V ) + 5T(x a ,n i ) (4) 

where 6T is much smaller that homogeneous component T . If the observer is at rest with respect 
to a certain coordinate system then taking into account that g a pu a u 13 = g 00 (u ) 2 = 1 we find 
that the photon frequency measured by this observer is equal to uj = Po/v^oo- ^ one S oes to 
the other coordinate system x a = x a + £° infinitesimally different from the old one, then the 
frequency of the same photon, measured by a different observer, who is at rest with respect to 
this new coordinate system, changes. From the transformation properties of the metric and the 
4-momentum one gets that 

co =>• u =p / v /^ = uj(1 + C' n*) (5) 

where I used the eq. p a p a = and kept only the first order terms in £ and metric perturbations; 
prime denotes the derivative with respect to time rj. Taking into account that the distribution 
function is a scalar, one easily finds that the small temperature fluctuations measured by an 
observer (at rest) in the new coordinate system are given by 

5T = 5T- + ToC'm (6) 

Hence, we see that only the monopole and dipole components depend on the particular coordinate 
system. The monopole component can always be removed by a redefinition of the background 
temperature and can not be measured locally. The dipole component depends on the motion of 
the observer with respect to the "new ether" defined by the background radiation and measuring 
it we can find how the Earth moves relative to CMB. Both of these components are not very 
interesting regarding the spectrum of the initial fluctuations. The higher multipoles depend 
neither on the particular observer or coordinate system we use to calculate them. Therefore I 
perform the calculations in conformal Newtonian coordinate system where these calculations look 
especially simple. 

Let us solve the Boltzmann's equation for the free propagating radiation in a flat universe 
with the metric 

ds 2 = a 2 {(1 + 2$) drf - (1 - 2$) 5 ik dx i dx k } (7) 
where $ 1 is the gravitational potential. Using the geodesic equations 

d\ P ' dX 2dx aPP ' { ) 

where A is an affine parameter, the Boltzmann's equation (j2J) takes the form 

df df 9$ df 

^ + „.( 1+ 2*)^ + 2p-^0. (9) 
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Taking into account that 

w = Po/v^oo = (! + $) p/a, (10) 

and using the Planck the ansatz ©, © one can easily get that in the lowest order in perturbations 
the Boltzmann's equation reduces to 

(T a)' = 0, (11) 

while the first order terms lead to 

d , d \ (ST \ 9$ 

In the most interesting case when the universe after recombination is dominated by dust, a 
nondecaying mode of the gravitational potential is constant and therefore the right hand side 
of the equation ()12j) vanishes. The operator in the left hand side is a total time derivative and 
therefore 

'ST 



— + $J = const, (13) 

along a null geodesies. The influence of the gravitational potential on the CMB fluctuations is 
known as Sachs- Wolfe effect. In the case when the gravitational potential is time dependent the 
combination (ST '/T + $) is not constant anymore. As it is clear from (fT2*|) its change is given 
by the integral from the partial time derivative of the potential along geodesies. This effect is 
usually called the integrated Sachs- Wolfe effect. If at late stages the universe is dominated by 
quintessence , the integrated SW effect can be essential, changing the resulting amplitudes of 
the fluctuations by 10 -j- 20 percent in big angular scales 6 > 1°. The accounting of this effect is 
rather obvious and therefore to avoid the overcomplication of the final formulae I consider only 
the case of the constant potential. 

As it follows from the geodesies equations the photons arriving at present time 770 to 
observer located at x l from the direction n l propagate along geodesies 

ar*(?7) ~ xl + n { (n - T] ) . (14) 

Therefore, from (|T3J) we get that ST jT in the direction n % on the sky is equal today to 

ST ST 

— (770,4,^) = ^(^ar^r),^) + [$(»&., ^(tfc.)) " $ {Vo, 4)] (^) 

where i] r is the recombination moment and x' l (r] r ) is given by (j!4jl . Since we live in a particular 
place we are only interested in n l -dependence of the temperature fluctuations. Therefore, the 
last term in (|15|) . contributing only to the monopole component, which is not measurable locally 
anyway, can be ignored. As we see the angular dependence of (ST/T) Q is determined by two 
factors: a) by the "initial value" of (ST/T) r in n— direction in a place from where the photons 
arrive and b) by the value of the gravitational potential $ in this place. The appropriate tem- 
perature fluctuations at the moment of recombination (ST/T) r can be easily expressed in terms 
of the gravitational potential and the fluctuations of the photon energy density 5 7 = Se^/e^ at 
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this time. With this purpose let us use the matching conditions for the hydrodynamical energy 
momentum tensor (EMT), which describes the radiation before decoupling, and the kinetic EMT 
(see, for instance, jH]) 

1 f rf^*P a PP 



Tp = —1= / dV*Hn ( 16 ) 

V-9 J P 

characterizing the gas of the free photons after decoupling. Substituting the expression for the 
metric into (JTfij) and using for the distribution function ansatz Q we get (up to the linear in 
perturbations terms): 

T ° * J d*pf (|) Po * T 4 J (l + 4^ f(y) y*dyd% (17) 

where y = uj/T and we have expressed po an d P through u. The integral over y from the Planckian 
function / can be explicitly calculated and give the numerical factor, which, being combined with 
47tTq, is equal to the energy density of the unperturbed radiation. Right before recombination the 
appropriate component of hydrodynamical EMT for the radiation is equal to Tq = e 7 (1 + 5 7 ) . 
This component doesn't jump at the moment when the universe becomes transparent and hence 

&y = 4 / — • (18) 

J T An v ; 

Similar by, one can derive from (fTBj) that for the kinetic EMT 

li -4* (19) 

On the other hand as it follows from the conservation law for the coupled photon-baryon plasma 
(jlH2j) (see appendix B) the appropriate divergence for the hydrodynamical components of Tq can 
be expressed in terms of 5 7 and $; hence 

where I have assumed that at recombination the cold matter dominates and therefore neglected 
the time derivative of the potential: $' (r] r ) = 0. Going to the Fourier space we easily infer that 



6 -p) k (Vr) = \(s k + ^ 2 (k m n m )S' k ) (21) 



satisfies both matching conditions (fTHj) and ([20)1 . Here and later on we skip the index 7 assuming 
that S always denotes the fluctuations of the radiation energy density. Taking into account the 
initial conditions (|2*T|) and skipping the monopole term in (|15|) . we obtain the following expression 
for the temperature fluctuations in the direction n = (n 1 ,n 2 ,n 3 ) at location xo = (a; 1 , a; 2 , x 3 ) 
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where k = |k| , k • n = k m n m and k • x = k n x^. Since r] r /i] < 1/30 we can neglect here r] r 
compared to r] . It is clear that the first term under the integral represents the combined result 
from the initial inhomogeneities in the radiation energy density and Sachs- Wolfe effect, while 
the second term is related to the velocities of the baryon-radiation plasma at recombination and 
therefore, is called Doppler contribution to the fluctuations. 



3 Correlation function and multipoles 

In the experiments one usually measures the temperature difference of the photons received by two 
antennae separated by a given angle 9 and this squared difference is averaged over the substantial 
part of the sky. The obtained quantity can be expressed in terms of the correlation function 

C(9) = ^(n^(n 2 )J) (23) 

where the brackets () denote the averaging over all ri! and n 2 , satisfying the condition rix • n 2 = 
cos(6 l ). Actually, 




T( ni )-r(n 2 " 2 

Tn 




2 (<7 (0) - <7 (0)) (24) 



On the other hand, for a given perturbation spectrum the correlation function C (8) can be 
easily expressed in terms of the expectation values of the Fourier components of the quantities 
characterizing these perturbations at the moment of recombination. 

The Universe is homogeneous and isotropic in big scales and therefore the averaging over 
the sky for a particular observer and a spatial averaging over the positions x should give for 
small angles (or big multipoles) nearly the same results. Therefore, the problem of averaging is 
finally reduced to the averaging of the products of Fourier components for the random Gaussian 
field. Substituting (J2~2~j) into (|2^j) and taking into account that, (c&k^V) = \&k\ 2 S (k + k') , after 
integrating over the angular part of k we obtain: 

h K d \ / , 6 k 35' k d 



C J [® k+ 4 Ak 2 d7 ll )y k+ A Ak 2 dr] 2 ) fc|n 1 77 1 -n 2 77 2 ] 2vr 2 ' ( " J ' 

where after differentiation with respect to rji and t? 2 we have to put rji = r] 2 = 1] . Now using the 
formula (see (10.1.45) in [0]): 

5i " ( f """ " n2 '| 2l) =f> + 1)A (k, h ) J, (tan) P, (cosO) (26) 
k |ni?7i - n 2 ?7 2 i 



1=0 
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where Pi (cos 0) and ji (kr]) are, respectively, the Legendre polynomials and spherical Bessel func- 
tions of order I, we can rewrite the expression for the correlation function in the following form 



oo 

C (0) = — ( 2l + 1) °i p i ( cos 0) (27) 

/=9 



Z=2 

where the monopole and dipole components (I = 0, 1) were excluded and 



c, = - 

7T 



®k C»7r) H j — J« ( fc7 7o) 



4 y J(V /u/ 4A; d(krio) 



2 

A; 2 d£;. (28) 



The coefficients C\ are directly related to the coefficients ai m in the expansion of 5T/T in terms 
of spherical harmonics, namely Ci = (\aim\ 2 ) , and therefore they characterize the contribution of 
the multipole component I to the correlation function. If 6* 1 the main contribution to C (0) 
give the multipoles with / ~ 1/0. 

The resulting spectrum of CMB-fluctuations depends of the various cosmological pa- 
rameters. First of all, these are the amplitude and the index of the primordial spectrum of 
inhomogeneities, generated by inflation. The rather generic prediction of inflation is that in the 
interesting for us scales 2 : \&lk 3 \ = Bk n ~ x , with 1 — n s ~ 0.03 0.08 0, 0- The amplitude B is 
not predicted and should be normalized to fit the observations. The other parameters on which 
the shape of the CMB-spectrum depends are the baryon density, characterized by f2&, the contri- 
bution of the clustered cold matter to the total energy density Q m (Q m = + ft c dm), the Hubble 
constant h 75 (I normalize it on 75 km /sec - Mpc) and the cosmological constant (quintessence) 
characterized by The present data are best fitted assuming that the universe is flat with 
Q to t = + — 1 and the total energy density is dominated by the dark cold matter and 
quintessence with only a few percent of baryons. Below I concentrate mostly on the models, 
which deviate from the "concordance model" not too much. Our purpose is to clarify how the 
variation of the parameters influences the observed CMB spectrum and to get an idea up to 
what extent the CMB determination of the cosmological parameters is robust. I will consider the 
different angular scales separately. 



4 Anisotropics in big angular scales 

The formula (|28j) was derived in the approximation of the instantaneous recombination. Because 
of causality this approximation is rather good when we consider big angular scales, where the 
CMB fluctuations are mainly determined by inhomogeneities exceeding the horizon scale at re- 
combination. Moreover, the perturbations spectrum in these scales is not much influenced by 

2 Most of the calculations will be done here for a flat spectrum (n = 1). The consideration can be easily 
generalized for an arbitrary spectral index n and in the case of small deviations from the flat spectrum the 
modification of the final results is obvious (it will be briefly discussed later). 
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the evolution. Hence the CMB fluctuations in big angular scales deliver us the undisturbed in- 
formation about primordial inhomogeneities, which are characterized by the amplitude of the 
primordial spectrum B and by the spectral index n s . The horizon at recombination is about 
the Hubble scale H~ l = 1.5t r , which in flat universe has the angular size 0.87° on today's sky. 
Therefore, the fluctuations which we will consider in this section refer to the angles 9 ^> 1° or, 
to the multipoles I < 1/9 H ~ 200. 

For the superhorizon adiabatic perturbations with kr\ r < 1 we have (see Appendix B): 

8 k (7j r ) ~ 5' k (7ft) ^ 0. (29) 

As it follows from ((221) their contribution to the temperature fluctuations is equal to 

ST 1 

— (r/ ,x ,n) ~ -$ (r/ r ,x - nr/ ) , (30) 

that is the observed fluctuations constitute one third of the gravitational potential in a place from 
where the photons arrived. Taking into account that after equality the potential in super curvature 
scales drops compared to its initial value $° by factor 9/10 3J, substituting (j2HJ) into (J2~%j) and 
calculating the integral with the help of the standard formula 

s Jl{s) - 2 7r P(^)r(/ + 2-f) {S1) 

($°) 2 A; 3 = B) we obtain well known result: 

9B 

(Z(Z + l)Ci), <30 = — = canst. (32) 



for the flat initial spectrum 



Deriving this formula I used in the integrand the flat spectrum everywhere, assuming 
that the main contribution for small I comes from the scales exceeding the horizon, where the 
primordial spectrum is not modified by evolution. This is a rather good approximation for I up 
to 20 -T- 30. Nonetheless, when we are interested in the precise normalization we need to take into 
account the corrections coming from the modified spectrum of the perturbations at big k. This 
can be well traced only in numerical calculations. 

Unfortunately, the accuracy of the direct information about the statistical properties of 
the primordial perturbations spectrum gained from the measurements in big angular scales is 
restricted by the cosmic variance. In fact, within cosmic horizon there are only 21 + 1 samples of 
the statistical realization for every particular multipole component /. This leads to the minimal 
inevitable typical "statistical fluctuations" in Ci 

^^{2l + l)- l '\ (33) 
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Hence, the statistical properties of the spectrum in the scales corresponding to the multipole 
/ can be determined in observations only up to an inevitable "error" given by (|33|). For the 
quadrupole (I = 2) this "typical error" is about 50% and therefore it can not be used for the 
normalization of the spectrum. For I ~ 20 the error constitutes 15%. Therefore, if we want to get 
a better accuracy in determining the spectrum of primordial inhomogeneities we are forced to go 
to smaller angular scales, where the spectrum is distorted by evolution. On the one hand it is bad 
news, since we lose the "pristine information". However, on the other hand, the distortions of the 
spectrum depend on the other cosmological parameters involving them "directly in the game" and, 
therefore, allowing us to determine these parameters under condition of having precise enough 
measurements. 

On small angular scales we can not ignore anymore the effect of the delayed recombination 
and the obtained above formulae should be corrected. Therefore before I proceed with calculations 
of the fluctuations in small scales I will find how the basic formulae should be modified to account 
for the effect of delayed recombination. 

5 Delayed recombination and finite thickness effect 

The delayed (non-instantaneous) recombination is important because of two reasons. First of all, 
the finite duration of recombination makes the moment when a specific photon decouples to be 
not very definite. As a result the information about the place from where this photon arrives 
is "smeared out". This leads to a suppression of the CMB-fluctuations in small angular scales, 
known as finite thickness effect. The delayed recombination leads also to an extra dissipation of 
the inhomogeneities increasing the Silk damping scale and hence changing the conditions in the 
places where the photons decouple. First we consider the finite thickness effect. 

Let us consider a particular photon arriving to us from the direction n. With non- 
negligible probability this photon could decouple at any value of the redshift in the interval: 
1200 > z > 900 and propagate without further scatterings afterwards. If this happens at the 
moment t)l then the photon arrives to us from the place x (tjl) = Xo + n (rji — r)o) without further 
scatterings and brings the information about conditions in this particular place. Since we do not 
know exactly when and where the particular photon decouples, a set of the photons arriving 
from a definite direction brings us only "smeared" information about the conditions within the 
layer of width Ax ~ Atjl, where Atjl is the duration of recombination. It is clear that if the 
perturbation has a scale smaller than Atjl then as a result of this smearing the information about 
the structure of this perturbation will be lost and we expect that the contribution of these scales 
to the temperature fluctuations will be strongly suppressed. 

Let us calculate the probability that the photon was scattered last time within the time 
interval At^ at the moment of physical time ti (corresponding to the conformal time t]l) and 
then avoided further scatterings until present time t^. With this purpose we divide the time 
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interval to > t > t^ into N small pieces of duration At, so that, tj = + jAt and TV > j > 1. It 
is obvious that the required probability is 

At^ / At \ / At \ ( At \ 



where r (tj) = (orn t (tj) X (tj))~ is the mean free time due to the Thompson scattering at t 
and rit, X are, respectively, the total number density of all (bounded and free) electrons and the 
degree of ionization. Taking limit N — > oo (At — * 0) and going back from the physical time t to 
conformal time rj we obtain: 

dP (r} L ) = //' (r] L ) exp (-/i (r) L )) dr] L , (35) 

where prime, as usually, denotes the derivative with respect to conformal time and the optical 
depth 

f to dt f vo 

A* (Vl) = / 777 = / <y T n t X e a (rj) drj. (36) 



was introduced. Now, taking into account that in the formula (|22jl the recombination moment 
r] r should be replaced by r/ L weighted with the probability (JHH1) we conclude that this formula 
should be modified as: 

j r 4 4^ o vo ) m e » exp ( /i) dr?i (27r ) 3 / 2 (37) 



T 



I would like to stress that in distinction from (|22|) here one can not neglect T)l compared to r] 
anymore since when we integrate over t]l the appropriate argument of the exponent changes very 
much for k > r^ L 1 . 

It is easy to see that the visibility function //exp (— fi) vanishes at very small rji (because 
fj, 3> 1) and at big T)l (//' — > 0) and reaches the maximum at r\ r determined by the condition 

// = /i' 2 (38) 

Since in the case of non-instantaneous recombination the moment when the photons decouple 
from the matter become smeared over rather substantial time interval we reserve from now on 
the notation r] r for the time when the visibility function takes its maximum value. This maximum 
is located within thin layer 1200 > z > 900. During this short time interval the scale factor and 
the total number density n t do not change very substantially and therefore we neglect their time 
dependence, estimating the appropriate values at rj = rj r . On the contrary, the ionization degree 
X changes by few orders of magnitude. Taking this into account we can rewrite the condition 
(JSH1) as: 

X' r ~ -{a T n t a) r X 2 r (39) 

where index r means that the appropriate quantities are estimated at rj r . At 1200 > z > 900 the 
ionization degree X is well described by the formula (|115|) in Appendix A. The change of X is 
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mainly due to the exponential factor there; hence 



X^- 1MX z Wi HX (40) 



where Ti = (of /a) . Substituting this relation in (|3*9J) we get 



X r ~ H r n (a T n t a) r 1 (41) 

where k = 14400/z r . Together with f)115[) this equation determines when the visibility function 
takes its maximum value. It is easy to see that this happens in the " middle" of the recombination 
layer at z r ~ 1050 irrespective of the values of the cosmological parameters. At this time the 
ionization degree X r is still k ~ 13.7 times bigger than the ionization degree at the moment of 
decoupling determined by condition t ~ r 7 (see (|119)) ). Near its maximum the visibility function 
can be well approximated by the Gaussian one: 

//exp (-/x) oc exp (ji - h\rf)' r (r} L - r] r ) 2 \ (42) 
Calculating the derivatives with the help of (jlDj) , (f4*T]) we obtain 

//exp (-„) - ^exp (nUnt (f r - A") (43) 



where the pre-exponential factor was taken to satisfy the normalization condition J // exp (— /i) dr]L 
1. 

We can use this formula to perform the explicit integration over t]l in (jHTj) . The grav- 
itational potential and the slowly varying contribution to 5 1 practically do not change during 
recombination. Therefore, they can be approximated by their values at r] r . The only term inside 
the curly brackets in ()37|) which could incur a very substantial change is the Silk damping scale. 
Keeping in mind that the main contribution to the integral comes from the region near r] r we 
estimate this scale also at r) r . Of course this is a rather rough estimate which nevertheless repro- 
duces the results of the numerics with rather good accuracy . Thus, ignoring r^— dependence of 
the expression in curly brackets in (j3T|) and taking its value at r] r , after substitution (f4l?|) in (jS7j) 
and integration over t]l we obtain 



T 

where 



d ' 1 ^ 4- - - ——\ PYD (- (nhn ?) ptt-Cxo+nfrr-m)) (A A) 



a = (45) 

VQ(KH V ) r 

In deriving I replaced (k ■ n) 2 by k 2 /3, keeping in mind the isotropy of the perturbations. 
Note that now we can neglect in the exponent r\ r compared to 770 . To find how a depends 
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on the cosmological parameters we have to calculate (Tirj) . At recombination and before the 
cosmological term can be ignored and the behavior of the scale factor is well described by solution 
for the matter-radiation universe (see, for instance, [3]), 

° ( '> =-((£)' + 2 (£))■ (46) 

Note, that 77* is a bit different from the equality time r] eq when the energy densities of radiation and 
matter are exactly equal. The relation r] eq = (a/2 — l) 77* follows from the equation a (r] eg ) = a m . 
Hence 

(*!>, = (47) 

2 + IVr/V*) 

where {f] r /f]*) can be expressed through the ratio of redshifts at equality and recombination if we 
use obvious relation 2 

+ 2 (V) ~ £f£ (48) 

J}* J \V*J Z r 

Substituting (|T7jl in (J43|) and taking into account that k ~ 13.7 we obtain: 

-1/2 N 

n — 1. 19 x 10- 2 [ 1+ (1 + ^ ) ) ( 19 ) 



The exact value of z eq depends on the cold matter contribution to the total energy density and 
from the number of the ultrarelativistic species. Assuming three types of neutrino z r /z eq can be 
estimated as 

^ ~ 12.8 (tt m h 2 5 ) (50) 

The value of a depends on the amount of the cold matter not very sensitively; if fl m h 2 5 — 0.3 
then a ~ 2.2 x 10~ 2 , while for tt m h 2 75 ~ 1, a ~ 1.9 X 10~ 2 . 

Now let us find how non-instantaneous recombination influences the Silk dissipation scale. 
As we mentioned above, at 77 = 7] r the ionization degree is k ~ 13.7 times bigger than at the 
decoupling and the mean free path of the photon is correspondingly smaller than the horizon 
scale. Therefore we can try to use the result (|141|) of Appendix B, obtained in imperfect fluid 
approximation, to estimate the corrections to the formula ()142j) due to noninstantaneous recom- 
bination. Using the approximate formula ()115j) which is valid when the ionization drops below 
unity we obtain 



(k DV ); 2 ~ 0.36 (n m h 2 75 ) 1/2 (fi 6 /4) - 1 z^' 2 + ^c 2 s a 2 (51) 

The first term here is the same as dissipation scale ()142j) derived for the case of instantaneous 
recombination. It accounts the dissipation until the moment when recombination begins. The 
second term is due to an extra dissipation which happens in the process of recombination. Note 
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that the second term in (|51|) corresponds to the scale which at t] r smaller that the mean free path 
r 7 and barely can be trusted literally. However, within the time interval Arj ~ T] r a when the 
visibility function is different from zero the free propagating photons have enough time only to go 
at the (comoving) distance A ~ rj r a which roughly corresponds to the second term in (j51|) . Hence, 
although the imperfect fluid interpretation of the second term becomes questionable, it can be 
nevertheless used to make an estimate of the damping scale. For the realistic values of the dark 
matter and baryon densities, Q m hj 5 ~ 0.3 and Qbhj 5 — 0.04, this term is nearly twice bigger than 
the first term; hence the extra Silk dissipation due to delayed recombination is rather important. 
At very low baryon density the first term in (joTj) dominates and most of the dissipation happens 
before ionization significantly drops. 

Thus, we found that the delayed recombination can be taken into account in a simple 
way. First, there occurs the extra dissipation of the perturbations and the dissipation scale can 
increase in few times compared to the case of instantaneous recombination. Second, it leads to an 
uncertainty when the photons decouple from the matter and as a result to an extra suppression 
of the CMB-fluctuations in small angular scales. Although both effects are interconnected they 
have different nature and should not be confused. 

The formulae derived in the approximation of instantaneous recombination are modified 
in an obvious way. Namely, the formula should be replaced by (JUJ). Repeating the steps 
which lead to the key formula (|28|) we conclude that the expression under the integral there 
should be just multiplied by a general factor exp (—2 (crkri r ) 2 ) . 

6 Small angular scales 

At big I, corresponding to small angular scales, the main contribution to C\ give the perturba- 
tions which being placed at recombination have an angular size 9 ~ 1/7 on today's sky . The 
multipole moment I ~ 200 corresponds to the sound horizon scale at recombination. Hence the 
perturbations responsible for the fluctuations with I > 100 4- 200 should have the wavenumbers 
k > rj^ 1 , that is, they enter horizon before recombination. These perturbations evolve in a rather 
complicated way and the primordial spectrum is strongly modified at k > i]~ x . In realistic models, 
the transfer functions relating the initial spectrum of gravitational potential $° with the result- 
ing spectra for $ and 5 1 at r\ r can be analytically derived only in two limiting ) for the 
perturbations which entered horizon well before equality and b) much later after equality (when 
the gravitational field of radiation can be ignored). In the limit of very big k the result is given by 
the formulae (|152|) . (J153|) . while for very small k by (|143|) (see Appendix B). Unfortunately, for 
the realistic values of the cosmological parameters none of these results can be directly used to 
calculate the CMB fluctuations in the most interesting region of first few acoustic peaks. Actually, 
the derived shortwave asymptotic is applicable only for those perturbations which have chance 
for at least one oscillation before equality (kr] eq > 2y / 3vr ~ 10). At the same time the longwave 
asymptotic ()143|) can be literally applied only to the perturbations which entered the horizon 
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when the radiation was already negligible compared to the matter. If Q m hj 5 ~ 0.3 then as it 
follows from (j3D)) z r /z eg ~ 4, and the radiation still constitute about 20% of the energy density at 
the recombination time. Hence, the formula (j!43|) is not trustable for those perturbations which 
enter the horizon in between equality and recombination and responsible for the fluctuations in 
the region of first few acoustic peaks. 



6.1 Transfer functions 

To describe the perturbations in these intermediate region we have to modify the derived formulae. 
Taking into account the time behavior of the asymptotic WKB-solutions of Appendix B we 
conclude that at the moment of recombination: 



8k 

+ — 

k ^ 4 



T p T oV /^cos (k jT c s d V j e- {k ' k ^ 2 



K (52) 



and respectively 



5' k ~ -AT kcl /2 sin [k c s dTi^J e ~ (k/kD)2 $ k (53) 

where the transfer functions T p and T should depend on the wavenumber k, equality time rj eq 
and baryon density Qb- To simplify the consideration we will restrict ourselves by the case when 
the baryon density is small compared to the total density of the cold matter, that is, Qb <C Q m . 
This will allow us to neglect the baryon contribution to the gravitational potential compared to 
the contribution of the cold dark matter, which interacts with the radiation only gravitationally. 
However even in this case the baryons influence the speed of sound and we have to take this into 
account. This is the situation for the concordance model and one can use analytical results, which 
I will derived below, only to study the dependence of the fluctuations on the values of the major 
cosmological parameters within some "window" around this model. If the contribution of the 
baryons to the gravitational potential is negligible the transfer functions T p and T Q depend only 
on k and t] eq , which on dimensional grounds can enter T p and T Q only in combination kr] eq . Their 
asymptotic behavior can be easily inferred from (|143|) . (|153|) . For the longwave perturbations 
with kr] eq <C 1, 

T p - — ; T -> — • 3- 3 / 4 ~ 0.4, (54) 
p 10' 10 ' V ' 

while in the shortwave limit for kr\ eq ^> 1 

InJOmgJ . r .^3* (55) 



(0.27kr] 



cqj 



where the factor 10/9 accounts for the change of the gravitational potential for superhorizon 
perturbations after matter-radiation equality. Unfortunately, in the most interesting intermediate 
range of scales 1 < ki] eq < 10 which is responsible for the fluctuations in the region of first 



15 



few acoustic peaks, the transfer functions can be calculated only numerically. In the interval: 
1 < ki] eq < 10, one can approximate T p with good accuracy bv[TU] 

T p ~ 0.25 In ] , (56) 



kr), 



eq 



and, respectively 3 , 

T ~ 0.361n(5.6fcr/ e ,). (57) 

The transfer functions are monotonic; as /c?7 e( j increases the function T p decreases and approaches 
zero, while T increases and reaches its asymptotic value T ~ 1.97. For perturbations which 
enter horizon well before equality, the function T Q is about five times bigger than for the per- 
turbations which cross the horizon late after equality. The physical origin of this difference is 
rather transparent. Before equality the gravitational field of the radiation can not be neglected. 
Therefore when perturbation enters horizon the gravity field of the radiation extra boosts the 
generated sound wave and its amplitude will be five times bigger than the amplitude in the case 
when the radiation can be neglected. 



6.2 Calculating the spectrum 



To calculate C\ we should substitute (JH2*|) . (JHH)) into formula (J2BI, which should be appropriately 
corrected for the finite thickness effect. However, the obtained integrals are not very transparent 
and before we proceed with their calculation, it makes sense to simplify these integrals using 
the advantage of considering I 3> 1. With this purpose we first get rid of the derivatives of the 
spherical Bessel function in (|28j) . Using the Bessel function equation one can easily verify that 



3? iz) 



1(1 + 1) 



J? (z) + 



(zjf (*))' 
2z 



(5* 



where prime denotes the derivative with respect to the Bessel function argument. Substituting 
this into (|28|) and integrating by parts we get 



k 2 + 



9\5'\ 
~L6 




l + 0)e 



-2{akt] r y 



jf(k Vo )dk (59) 



where by O I denoted the corrections of order r] r /r]o and (krjo) , which were estimated taking 
into account the general structure of the expressions in (|52"|). The corrections rj r /rj can be 

neglected compared to unity since rj r /rj < z r ~ 1/30. At big I only those k give a substantial 

3 I would like to thank A. Makarov for numerical calculations of T function in the limit of vanishing baryon 
density. 
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contribution to the integral for which kr]Q > I. Actually, as I — > oo we can use the following 
approximation for the Bessel functions 

31 (Z) ~* { z- 1 ' 2 (z 2 - z/ 2 )" 1/4 cos - v arccos [yjz) - vr/4) , z > u, ^ 

where vjz 7^ 1 is held fixed and v = I + 1/2; hence the correction l/krj ~ 1/Z <C 1 can also be 
skipped. 

Now I will use (IfjQjl in the integrand of (|59jl . Keeping in mind that the argument of jf (krjo) 
changes with k much faster that the argument of the oscillating part of the WKB solutions for 
5k let us replace the cosine squared, coming from (J5U|) . by its average value 1/2. The result reads 



c ,4f |4 y |V ^ ^ (61) 

where using the advantage of considering only big multipoles I replaced I + 1 with /. This result 
was first derived in [7j. 



Let us consider the flat initial spectrum: 



B. Substituting (JS2J), 452) into (JET) 
ementary calculations we arrive to the 



and changing the integration variable to x = fcr/o/Z after e 
following result: 

l 2 Q ~-(0 + N) (62) 

7T 

where keeping in mind /—dependence of l 2 Ci I have written it as a sum of different terms. Namely, 

= 0, + 2 , (63) 
is the oscillating contribution to the spectrum given by two terms with twice different periods: 

v V Ji x 2 V^T 

and 

^ {l-9c 2 )x 2 + 9c 2 _ nnM 



2 = f / ? T ^s^_r^ e - ms y x > cos {2lQx) ^ (65) 

These terms modulate the spectrum, leading to the peaks and valleys. I have introduced here 
the ratio 

1 f Vr 

9 = — c s [T]) drj, (66) 
Vo Jo 

which determines the period of oscillations and location of the peaks. The scales If and lg 
characterizing the damping of the fluctuations because of the Silk dissipation and finite thickness 
effect are equal to: 

lj 2 = 2a 2 ^- lf = 2(a 2 + (k DVr )- 2 )(^\ (67) 
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where a is given in (J49j) . The analytical estimate for the Silk scale kr>T] r is not very accurate, 
however one still can use the estimate (|51|) for k^Vr- 

In turn, the nonoscillating contribution J c can be written as a sum of three integrals 

N = Ni + N 2 + N 3 , (68) 

where 

/ , \ 2 „oo rp2 e -(l/lf) 2 x 2 

N x = 1 - / p , dx (69) 

is proportional to the baryon density and vanishes in the absence of baryons when c 2 = 1/3. The 
other two integrals are: 

roc T 2 -(l/lgfx 2 

N2 = i 2r^-r dx > ( 7 °) 
1 J i iVx — 1 

and 



N z= 9 4f T 2 ^_±e-^* 2 dx. (71) 

Before we proceed further with the calculation of the integrals let us express the parame- 
ters entering (jo^j) . namely, c s , If, Is, Q and transfer functions T Q , T p through the basic cosmo logical 
parameters fib, ^m, h 75 and = 1 — Q m . 



6.3 Parameters 

The speed of sound c s at recombination depends only on the baryon density, which determines 
how it deviates from the speed of sound in purely ultrarelativistic medium. To characterize these 
deviations it is convenient instead of the baryon density to introduce the parameter £ defined as 

ts ^- 1 = i{Sr 17 ^' (72) 

Then c 2 s can be expressed through £ as 



s 3(1+0 

For the realistic value of the baryon density Qbhj 5 ~ 0.035 one gets £ ~ 0.6. 

The damping scales If, Is are given by (|67j). It is clear that to express them through 
the cosmological parameters we first have to calculate the ratio f] r /vo, which also depends on 
the cosmological term. To calculate this ratio let us consider an auxiliary moment of time 
?7o > Vx > Vr, when the radiation is already negligible and the cosmological term is still not 
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relevant for dynamics. Then to determine f] x /r] Q we can use the exact solution describing a flat 
universe filled by the matter and cosmological constant: 

/ 3 \ 2/3 

a (t) = a I sinh - H t j (73) 



As a result we obtain: 

(74) 



(sinh re) 



2/3 



with the upper limit of integration y = sinh 1 (Vt^/Vtm) 1 ^ 2 • Taking into account that fi^ = 1 — Q m 
one can use the following numerical fit for /a in ([74)1: 

J A ~ fi-°- 09 ( 75) 

which approximates the exact result with the accuracy better that 1% everywhere within the 
interval 0.1 < fl m < 1. 



The ratio r] x /r) r can be calculated with the help of ({Hi}) and is equal to 



" / l + 2^ ^ =(^-) 1/2 ((i + ^l] 1/2 (7(3) 



where we used the equation (|48jl to express ?y*/?7 r in terms of z eq /z r . Combining this formula with 
(174*)) we obtain 

* 1 (l l + ^Y-{-T\h (77) 



V \ \ ^eq / \ ^eq 

Substituting this together with the expression (jjQ)) for cr into (|*T7)) one gets 

1/2 

r \ 

Z, 



•'eg 



J, ~ 1530 ( 1 + 5- ) II 1 (78) 



where the ratio of the redshifts at recombination and equality for three neutrino types (see (|50*i ) 
is equal to 

7.8 x 10~ 2 (ttrnh^y 1 (79) 



7* 



^eq 



The scale // characterizes the damping of CMB-fluctuations because of finite thickness effect. It 
depends on both cosmological term and fl m h% 5 not very sensitively; for instance, if Q m hj 5 — 0.3 
and fi A /iy 5 ~ 0.7 we have If ~ 1580, while for £l m hj 5 ~ 1 and ^^75 ^ one gets // ~ 1600. 
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The scale Is describing the combined effect from the finite thickness and Silk damping 
can be calculated similar by. Using the estimate (|51|) for Silk dissipation scale one can easily find 
that 



/ 



0.71 , 



1 + 0.56£ 



+ 



0.8 



(^m^75 



,V2 



\ 



-1/2 



£(1 + ( 1 + {1+Zeq/Zr y 



1/2 



JO) 



This formula is not as reliable as the estimate for If since first we neglected the contribu- 
tion of the heat conductivity to Silk dissipation scale and second we calculated it using imperfect 
fluid approximation which surely breaks down when the visibility function reaches its maximum. 
Nevertheless it is still trustable within 10% accuracy and an exact result is a bit smaller than 
given by (|80|) . In distinction from If the damping scale 1$ depends not only on the matter density 
and cosmological term but also on the baryon density, characterized by £. However, this depen- 
dence is very strong only for £ <C 1 when the second term inside the bracket in (JHUj) dominates. 
For £ = 0.6 we get Is — 1100 if VL m hj b ~ 0.3 and Is — 980 for £l m h 75 ~ 1. The dissipation scale in 
the universe with more cold matter is bigger (correspondingly Is is smaller) because in this case 
the recombination happens at later cosmic time t r and hence the perturbations get an extra time 
to be washed out. 

The parameter p, which determines the location of the peaks, can be easily calculated if 
one substitutes the speed of sound 

where a{rj) is given by (|46|). into (JBTIj) and performs an explicit integration there. The result is 

i a , ( ^{i + z r / Zeq )l+ yTT+O 

Q — ~i= m = 

\ 1 + V£ (Zr/Zeq) 

It is clear that g depends on both baryon and matter densities. However, it is not very transparent 
how g behaves when we change these parameters. Therefore it is worthwhile to find a simple 
numerical fit for (|82|) . which would reproduce the parameter dependence of g within reasonable 
range of change of £ and Vt m h% b . The fit 4 

g ~ 0.014 (1 + 0.13O" 1 (tomffis) V4 J a (83) 

reproduces the exact result ([820 with the accuracy about 5 4- 7% or better everywhere in the 
region < £ < 5, 0.1 < Q m h 75 < 1 where the function g itself changes in about three times. 
Combining this with the numerical fit for Ja in (j75|l we have 

g ~ 0.014 (1 + 0. 130 _1 (n m h 3 7 i)°- W (84) 
4 I am very grateful to P. Steinhardt for helping me to check numerically the accuracy of this fit and the fits 
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The transfer functions T p , T a depend only on kr\ eq and can be expressed as the functions 
of variable x = krjo / 1 : 

kr]eq = !hii x ~ 0.72 (n m h 2 75 ) ~ 1/2 I A l 200 x (85) 

^0 

where Z200 = 1/200- As we will see the contributions to the integrals defining the fluctuations in 
the region of the first few acoustic peaks comes from O (1) > x > 1. Therefore for 200 < I < 1000 
the transfer functions in the relevant range of krj eq can be approximated by (|56p. ()57j) ; hence 



where 



and, respectively, 



T p (x) = 0.74- 0.25 (P + lnx) 
3 {I, tt m , h 75 ) = In [ ^ ) /20 ° 2 ) , 
T (x) = 0.5 + 0.36 (P + lnx) 



(86) 
(87) 
(88) 



6.4 Calculating the spectrum (continuation) 

Now I will proceed with the calculations of the fluctuations. The main contribution to the 
integrals from the oscillating functions (|64j) . (J65j) gives the vicinity of the singular point x = 1. 
These integrals have the form 

/(x)c ° s(ax) rfx (89) 



1 y/x — 1 

and after making substitution x = y 2 + l can be calculated using stationary (saddle) point method. 
The result is 

f (x) cos (ax) , f (1) Pk f 7T 1 D \ , , 

J dx& M ,, t -cos a+- + -arcsm _ , (90) 



yfr=l (1 + B 2 ) 1/4 V a V 4 2 VI + D 2 

where D = (c?ln//aiia:) a . =1 . For big a we can put D m and the above formula simplifies to 



f OC f(x)cOs(ax) J pK ( 7T 



(91) 



Using (}9*Tj) to calculate the integrals in (jHU), we obtain: 

O ~ (A cos (/^ + tt/4) + A cos (2Zg + tt/4)) (92) 
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where the coefficients 

are slowly varying functions of Z. They also depend on the basic cosmo logical parameters and the 
spectrum of the fluctuations at I > 200 is rather sensitive to the variation of these parameters. 
It is worth to mention that in this approximation the contribution of the Doppler term to the 
oscillating part of the spectrum drops out. One can check that actually this contribution at 
/ > 200 do not exceed few percent of the total amplitude. If fib <C fl m the transfer functions for 
the most interesting range 200 < I < 1000 can be approximated by (jHSj) and (|S5jl. In this case 
we have 



((P-0.78) 2 -4.3)£ r^/- 2 r a W „ (0.5 + 0.36P) 2 
Ai ^ 0.1^ '—j-, ^eUt'a -V J 1 J, ^ 2 ~0.14^ TW ^- 

(i + 1/4 (i + 1/2 



(94) 



where P is given by IJB7jl 

Substituting (JSfij) in the expression ()69|) for non-oscillating contribution Ni we get 

Ni ~£ 2 [(0.74 - 0.25P) 2 J - (0.37- 0. 125P)/i + (0.25) 2 / 2 ] (95) 

where the integrals 

can be calculated in terms of the hypergeometric functions. However the obtained expressions 
are not very transparent and therefore it makes sense to find the numerical fits for them. The 
final result is 



„ (P- 0.22 (1 l f a3 -2.6) /,,, y» 

Ni ~ 0.063£ 2 ^ ra — '-e-W 1 ') . 97 

1 + 0.65(Z/Z/) 1,4 



Similar by, we obtain 



0.037 (P - 0.22 (Z// 5 )°- 3 + 1.7)' _ 



1/2 1 + 0.65 (Z/Z^ 1 ' 4 



e-(^) 2 (98) 



The Doppler contribution to nonoscillating part of the spectrum is comparable to N 2 and is equal 
to 

N3 . (99) 



1 + £)' /2 1 + 2 (Z/Z 



The numerical fits ()97 )) -([99 )) reproduce the exact result in the most interesting range of multipoles 
with a few percent accuracy for a wide range of cosmological parameters. The extra dependence 
on I /Is and I /If is due to the fact that the exponent in the integrals from nonoscillating functions 
can not be just simply estimated at x — 1. When the expression under the integral is monotonic 
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function the substantial contribution to it comes not only from the vicinity of x = 1 but also from 
x ~ 0(1). The nonoscillating contribution of the Doppler term given by N 3 is rather essential 
and can not be ignored. 

It is convenient to normalize 1(1 + 1) G\ for big / to the amplitude of fluctuations for small 
I, given by ()32j). so that finally we obtain 



i(i + i)Q ioo 

(l(l + l)Q) l<30 9 1 



-N 1 + N 2 + N 3 ) . 



(100) 



where O, JVi, N 2 , N 3 are respectively given by (|H2J), (|9*7j). (|9*Hj). (}9*9"j) . In the case of the concordance 
model (Q m = 0.3, I^a = 0.7, Qb = 0.04 and = 70km/ sec ■ Mpc) the result is presented in 
Fig.l, where I have separately shown by the dashed and thin solid lines, respectively, the overall 
nonoscillating and oscillating contributions. The total resulting fluctuations are shown by the 
thick solid line. 
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About accuracy. Comparing (jlOOj) to CMBFAST runs 5 one can easily check that the 
analytical approximation works rather well reproducing the numerical results with good accuracy 
in a rather wide range of the cosmological parameters around concordance model. Namely, for 

5 I am very grateful to P. Steinhardt, S. Bashinsky and U.Seljak for performing numerous CMBFAST runs 
necessary for this work. 
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Vt m = 0.3 the agreement is still very good up to f2& ~ 0.08, when the baryons constitute already 
about 30% percent of the total cold matter density. At higher Q b the contribution of the baryons 
to the gravitational potential which we neglected becomes very essential and one can not use 
anymore the analytical formula (jlOOj) . This formula was derived under assumption <C Q m and 
is not trustable when the baryons constitute a very substantial fraction of the total amount of 
cold matter. It is also worth to mention that at high VL m h^ 5 the expected accuracy in the region 
of the first peak is not as good as in the region of the second and third peaks . This is because the 
used approximations for transfer functions become not so accurate on the border of the interval 
corresponding to 1 < ki] eq < 10. In particular, if Q m h% 5 = 1 the main contribution to the first 
acoustic peak located at I ~ 200 give the perturbations with the wavenumbers kr\ eq ~ 0.7 (see 
(|85|)) where the approximations (J56|) and (J57|) are not very accurate. Hence, although for the 
model with Q m hj 5 = 1 and Qf, = 0.04 the analytical result is still in a fair agreement with the 
numerics, its accuracy in the region of the first peak is not as good as for concordance model. Also 
note that the peaks given by (jl00|) are shifted by about 10% compared to the numerical results. 
One of the reasons for that is that with the purpose to simplify the final expression we neglected 
in (J§Uj) an extra phase shift proportional to D. The other reason is that we underestimated the 
parameter g which was derived in the assumption of instantaneous recombination. In reality the 
recombination takes place within about quarter of the cosmological time and in the process of 
recombination the baryons decouple from the radiation. As a result the sound speed increases 
and q should be a bit bigger compared to (f%2~j). 

However the main value of the analytical result is not in its competitive accuracy with 
the numerics, but because it allows us to understand the main features of the CMB spectrum and 
study explicitly how they depend on the cosmological parameters. In turn it opens a possibility 
to understand the degeneracy of the spectrum with respect to the certain combinations of the 
cosmological parameters which could lead to a "cosmic confusion". 

7 Determining the cosmological parameters 

Let us discuss how the main features of the spectrum change when the cosmological parameters 
vary. These parameters are: the amplitude B and the slope n s of the primordial spectrum, the 
baryon density characterized by the total cold matter density Q m , the cosmological constant 
Q\ and the Hubble constant h 75 . The amplitude and the slope of the primordial spectrum can 
already be determined with a reasonably good accuracy when we consider only the measurements 
in big angular scales. From these observations it follows that the spectrum does not deviate too 
much from the scale invariant (n s = 1). The cosmic variance, which is important in the big scales 
which are not much disturbed by the transfer functions does not allow us to conclude anything 
about small deviations from scale invariant spectrum predicted by inflation on basis of only these 
observations. 

The results for the fluctuations in small angular scales were derived assuming a flat 
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universe, where Vt m + Vt\ = 1, and a scale invariant spectrum with n s = 1. How do they change 
when the spectrum deviates from the scale invariant will be pointed out below. First, I would 
like to concentrate on flat models with scale invariant spectrum (n s = 1) and find how the 
characteristic features of the CMB spectrum depend on the cosmological parameters Qb, &m and 
h 75 (the cosmological constant is fixed by the flatness condition to be = 1 — fl m ). 

The location of the peaks and the flatness of the universe. The most interesting feature 
of the spectrum is the presence of the peaks and valleys, the height and location of which very 
sensitively depend of the major cosmological parameters. At / > 1000 the fluctuations are 
strongly suppressed and therefore the most interesting part of the spectrum is those one where 
the first three peaks are located. These peaks arise as a result of superimposing of the oscillating 
contribution to the fluctuations O, given by (|92j) . on the "hill" iV (I) = Ni + N 2 + N s representing 
a nonoscillating part of the spectrum (see Fig.l). It is clear that the locations and the heights of 
the peaks depend not only on the oscillating part, but also on the shape of the "hill". Let us 
neglect for a moment the effect of the "hill shape". In this case the location of the peaks would 
be determined by the superposition of two cosines in (P?2j). If |Al| <C A 2 the peaks should be 
located at 



where n = 1, 2, 3... and g is given by (|HB*j) . The first term in ()92|) has twice bigger period than the 
second and its amplitude A\ is negative. Therefore it participates in the constructive interference 
for the odd peaks (n = 1,3, ...) and in destructive interference for the even peaks (n = 2,4, ...). 
Moreover, because of the shift of the arguments of two cosines, the maxima of these cosines do 
not coincide and, as a result, first and third peaks (for which the interference is constructive) 
should be located in between the appropriate maxima of these two cosines, that is, at 



where the symbol -f- denotes the appropriate interval. If |^4i| ^> A 2 the peaks move closer to 
the lower bounds of the intervals in ()102j) . In fact, the situation is more complicated because the 
nonoscillating contribution N is not constant but is represented by "hill" . As it is clear from Fig.l, 
this leads to the further shift of the peaks to the "top of the hill" . For instance, for concordance 
model first peak moves a bit to the right, while the third peak to the left. Substituting £ ~ 0.6 
and Q m hj 5 ~ 0.26 into ()82|) we find that for this model the first peak should be located at 
l\ ~ 225 -T- 265, that is, somewhere in between 225 and 265. For the third peak l 3 ~ 825 4- 865. 
Because of the reasons I mentioned above, this result should be corrected by about 10% shifting 
the first peaks to the left. 

In the region of the odd peaks one has destructive interference of the oscillating terms. 
The first term in (|H2*|) which takes the minimal (negative) value tries to annihilate these peaks. 
The second peak (if it exists), should be located at 




(101) 




(102) 




(103) 



25 



or in the concordance model at I2 — 525 4- 565. 

How sensitive is the peaks location to the variation of the cosmological parameters? 
According to (|82p (see also (|84|)) g changes when the baryon and cold matter densities vary and 
therefore (see (|1U2|) . (jl03j0 the peaks location should also depend on these parameters. The 
parameter g is not very sensitive function of Q m , /175 and £. Therefore, the location of the first 
peak in a flat universe is relatively stable when we vary these parameters. In particular, when 
the baryon density increases in two times (£ ~ 0.6 — ► £ ~ 1.2) the first peak moves to the right 
by All ~ +20 and the shift of the second and third peaks are, respectively, A/2 ~ +40 and 
Al 3 ~ +60. When determining the location of the peaks, the baryon density always enters in 
combination £ oc fi&Ai^ with the Hubble constant. The cold matter density comes together with 
h 75 as Q m hf^. The increase of the cold matter density has an effect opposite to the increase of 
the baryon density, namely, if for given £ ~ 0.6 the cold matter density increases twice (fl m hfl ~ 
0.3 — > Q m hj^ ~ 0.6), the first peak goes to the left by All ~ —20 and respectively Al 2 ~ —40 and 
A/3 ~ —60. Thus we see that even in a flat universe, one can shift the location of the first peak 
quite substantially (All ~ 40) increasing the baryon density twice and simultaneously decreasing 
the cold matter density by the same factor. 

Why in this case can we be sure that the first peak location is a good indicator of the 
universe curvature? Fortunately, if we will fix the height of the first peak, then its location 
becomes "stable" with respect to the admitted variations of the cosmological parameters. The 
height of the first peak sensitively depends of the cold matter and baryon density. Given the height 
of this peak we can still vary the baryon and cold matter densities together. However if the cold 
matter density would increase and we would like still to keep the height of the peak to be the 
same, we have simultaneously change the baryon density, namely, it also should increase. Since the 
change of the baryon and cold matter densities have opposite effects on the peak location, it will be 
shifted not very much if both of them will increase simultaneously. For instance, if both densities 
increase by a factor two around concordance model, one can expect that All ~ 0. This explains 
the stability of the location of the first peak for the acceptable range of change of the cosmological 
parameters in a flat universe. The obtained result on the location of the first Doppler peak and 
its relative stability to the variation of the unknown cosmological parameters is a fair agreement 
with numerical calculations. The stability of the first peak location makes it an irreplaceable 
indicator of the total energy density of the universe. Actually, the peak location is incomparably 
more sensitive to the total energy density (in the open universe without cosmological constant 

— 1/2 

li oc Q tot ). The present observations strongly favor a flat universe (fi toi = 1), as predicted by 
inflation. 

Height of the peaks and the baryon and cold matter densities. In concordance model the 
amplitude of the first acoustic peak is in about 7 4 8 times bigger than the amplitude of the 
fluctuations in big angular scales. Substituting l n , given by (jl02|) . (|1U3|) . into (|57|) and using the 
formula (J82)) for g we see that the factor I\ is cancelled in the expression for P and therefore the 
height of the peaks given by (|100|) estimated at l n can depend only on fl m hj 5 and flbh% 5 (or £). 
If, for fixed Q m hj 5 , one increases the baryon density, the height Hi (Q m h? 5 ,£) also increases. In 



26 



the concordance model the increase of the baryon density by factor two (from £ ~ 0.6 to £ ~ 1.2) 
leads to the increase of the amplitude Hi in 1,5 times. This increase in the amplitude is mostly 
due to two terms in (jlU0|) . Ni (proportional to £ 2 ) and O (since A\ oc £). In turn, the increase of 
the cold matter density (at fixed £) suppresses the height of the first peak Hi. It becomes clear 
why this happens if we note that for fixed I, the function P entering the formulae for fluctuations 
decreases when Q m h% 5 increases. As a result the overall amplitude of the first peak decreases 
(mainly because N 2 and N 3 contributions decrease when VL m ti^ b increases). Therefore the height 
of the first peak is degenerate with respect to a certain combination of the baryon and cold 
matter densities. In a certain range of parameters the increase of the height due to the baryon 
density can be compensated if we simultaneously increase the cold matter density. However, if the 
baryon density would be too high, the increase of the height of first peak could not be anymore 
compensated by increase in can not much exceed unity. (Moreover, for 

big Q m hj 5 the transfer functions responsible for Q m hj 5 — dependence of Hi reach their asymptotic 
values for those values of krj eq which mainly contribute to the fluctuations in the region of the first 
peak). Hence, just relying on the result about the height of the first peak one can safely conclude 
that the baryon density can not be more than 15 -j- 20% of the total critical density. 

The degeneracy in determining Q m hj 5 and £ parameters can be easily resolved if we 
consider the second peak, which results mostly from the destructive interference of the oscillating 
terms in (I92J) superimposed on the "hill" given by iV-contribution. In the concordance model 
this peak is strongly suppressed in O— contribution and partially recovered only in the resulting 
spectrum because of the N— contribution (as one can see in Fig.l the "hill" has a sufficiently 
steep decline in this region). The presence of the second peak sensitively depends on the ratio 
of the amplitudes A\ and A?.- Since the amplitude A\ of the first term in (|§2~j) . which tries to 
"kill" the peak is proportional to the baryon density £, while Ai slightly decreases when £ grows, 
one can expect that the presence of large amount of baryons should diminish and may be even 
completely remove the second peak. Actually the "O— contribution" to the peak disappears when 
the baryon density increases only twice compared to its value in concordance model. However, 
in the resulting spectrum this peak still survives. This is because the growing amount of baryons 
simultaneously amplifies N\ — contribution to nonoscillating part of the spectrum and in turn this 
significantly steepens the "hill" in the region where the second peak is located. The analytical 
formulae become inapplicable at very high baryon densities. However the numerical calculations 
show that for Q m hj 5 ~ 0.26 the second peak is still present and has nearly the same amplitude 
as the third peak even if baryons constitute about 70% of all cold matter. Hence the presence 
of the second peak can not alone be considered as the indication of the low baryon density. 
Nevertheless, in combination with the observed height of the first peak the second peak is a very 
sensitive indicator not only for the baryon density, but also for total cold matter density. Given 
the height of the first peak, we can still vary the baryon and cold matter densities increasing or 
decreasing them simultaneously, since they "act in opposite directions" . However they influence 
differently the second peak. Namely, the simultaneous increase of the baryon and cold matter 
densities tries to "annihilate this peak". Actually, the amplitude of the second peak depends 
on the amplitudes Ai (£, fl m h% 5 ) and A2 ((,,Q m hj 5 ) in superposition of two cosines in (J02J)- The 
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increase of the baryon density tends to "kill" this second acoustic peak. The increase of the cold 
matter density at fixed £ has a similar effect. This is because A2 oc (T^) x=1 decreases faster than 
A\ oc (T T p ) x=1 when Q m hj 5 increases. At big Q m h% 5 the term which "kills the peak" dominates. 
Hence the height of the second peak depends simultaneously on the baryon and total cold matter 
densities and is very sensitive to the independent variation of both of them. Fixing the relation 
between £ and Q m h% 5 from the height of the first peak we can find the particular values of these 
parameters measuring the height of the second acoustic peak. For instance, if Q m hj 5 = 1, then 
£ should be about unity if one want to get the height of the first peak to be in agreement with 
observations. In this case second peak completely disappears. Hence the experimental detection of 
the second peak proves that the total density of the cold matter is smaller that the critical one and 
the baryon density is smaller than 6 4- 8%. This is in an excellent agreement with nucleosynthesis 
bounds. Shortly both of the results could be formulated as "too much baryons would destroy 
all deuterium and kill the second acoustic peak". In combination with the location and height 
of the first peak the presence of the second peak is also a strong independent indicator of the 
dark energy in the universe. In fact, from the location of the first peak it follows that the total 
density in the universe is critical and the presence of the second peak means that the cold matter 
can constitute only the fraction of it. 

Since the heights and locations of the peaks depend on the different combination of Q m 
and /i75 this allows us to resolve the degeneracy in determining the Hubble constant. As we 
have seen for a given £l},h% 5 the location of the peaks depends on Q m hf^, while their heights is 
determined by Q m h% 5 . Therefore keeping VL h h^ and Q m h% 5 to be fixed by the heights of the peaks 
we can still vary the Hubble parameter /175 shifting the position of the peaks. As it follows from 
(jl(J2|) . (jBU), for the given Vt b h^ ~ 0.04 and fi m /if 5 ~ 0.3 the increase of the Hubble constant by 
20% (say from 70 km/ sec ■ Mpc to 85 km/ sec ■ Mpc) moves the peaks to the left by 3%, that is, 
All — 7 and Al 2 — 15. Hence if we want to get an accurate determination of the Hubble constant 
from CMB spectrum alone we have to know the location of the peaks with very high accuracy. 
If the location of the peaks will be determined with 1% accuracy then the expected accuracy of 
the Hubble constant will be about 7%. 

Up to now we were assuming that the primordial spectrum of the inhomogeneities is 
scale invariant, that is the spectral index is n s = 1. The inflation predicts that there should be 
deviations from the scale invariant spectrum and we expect that n s ~ 0.92 -f- 0.97. The above 
derivation for the CMB fluctuations can be easily modified to account for these deviations. 

If n s 7^ 1 the obtained amplitudes of the fluctuations at given / should be just multiplied 
by the factor proportional to l l ~ n . To resolve the degeneracy in determining the cosmological 
parameters in this case the heights and location of the first two peaks are not sufficient. Actually 
for a given n one can always find the combination of the flbhy 5 and Q m hj 5 — parameters to fit 
the heights of the first two peaks. The location of these peaks is also not very sensitive to the 
deviations of the spectral index from unity. Therefore one needs extra information. With this 
purpose we can use for instance the height of the third acoustic peak. As one can check the 
height of this peak is not so sensitive to flbh% 5 and Vt m h^ as for the first two peaks. Fixing these 



28 



parameters and varying the spectral index n s for a given unchanged height of the first peak (this 
can always be done if together with n s we vary the amplitude of the spectrum B) we find that 
the relative height of the third peak changes as 



For instance, if n s ~ 0.95 the height of the third peak increases by about 5% compared to the 
case of n s = 1. From this estimate one can get a rough idea about necessary accuracy of the 
measurements to find the expected deviations from the scale invariant spectrum. 

8 Conclusions 



In this paper I have shown that if we assume that the main ingredients of the cosmological model 
are known, then we can completely resolve the degeneracy and determine the main cosmological 
parameters from the CMB spectrum. For that we just need to know the main features of the 
spectrum, namely, the heights and location of the peaks. Of course, the accuracy of the deter- 
mination is different for different parameters and seems to be the worst for the Hubble constant. 
The information we gain in the observations exceeds the discussed features of the spectrum. 
Namely, one measures also the entire shape of the spectrum, which, of course, also depends on 
the cosmological parameters. The necessity to fit this shape restricts the possible values of the 
parameters even in the case when we have the measurements only in the region of the first peak. 
This shape (as well as the heights and location of the peaks) also depends on the dissipation 
scales If and Is, which in turn slightly depends on the cosmological parameters. For the concor- 
dance model Is ~ 1000, and it is clear that the dissipation does not influence very much the first 
peak and becomes very essential in the region of the second peak and at high /. In particular, at 
I > 1000 this effect entirely dominates, leading to the exponential falloff of the spectrum at very 
high multipoles. This falloff is very sensitive to the parameters and, being measured, can give us 
extra information about them. The measurements of the polarization provides additional valu- 
able information about the cosmological parameters. When we vary the parameters the detailed 
behavior of the spectrum is, of course, more complicated than I described above (I also neglected 
here the primordial gravity waves which can give rather substantial contribution at / < 30). 
However, the above consideration correctly reflects the main features of this behavior and gives 
the physical understanding why the CMB spectrum so sensitively depends on the cosmological 
parameters. 
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A Hydrogen recombination 

The equilibrium description of recombination by Saha's formula fails nearly immediately after the 
beginning of recombination when only a few percent of hydrogen becomes neutral. Therefore one 
has to use the kinetic approach to describe the noninstantaneous (delayed) recombination [TT] . 

The direct recombination to the ground state with the emission of one energetic pho- 
ton is not very efficient. The emitted photon has enough energy to immediately ionize the first 
neutral hydrogen atom it meets. One can easily check that the two competing processes, direct 
recombination to the ground state and ionization, occur with a very high rate leaving no net con- 
tribution. More efficient is the cascade recombination when the neutral hydrogen is first formed 
in the excited state and then goes to the ground state. However, even in the cascade recombina- 
tion at least one very energetic photon is emitted. Its energy corresponds to the energy difference 
between 2p— and Is— states. This Lyman-alpha photon (L a ) has the energy 3.B#/4 =^> 117000°/^ 
and a rather big resonance absorption cross-section which at the recombination temperature is 
about a a ~ 1CT 17 -T- 1CT 16 cm 2 . The L a — photons are reabsorbed in r a ~ (cTann) 1 ~ 10 3 -j- 10 4 
sec after emission. This time has to be compared to the cosmological time. During the matter 
dominated epoch the cosmological time can be easily expressed through the temperature if we 
just equate the energy density of the cold particles to the critical energy density e cr = 1/ (67rt 2 ) 
and note that T = T 70 (1 + z) ; hence 



where is the Hubble constant normalized on 75 km/ sec ■ Mpc, T 7 o ~ 2.72 K and Q m is the 
contribution of cold matter to the critical density. At the moment of recombination r a <C t c ~ 10 13 
sec and the L a — quanta are not significantly redshifted before they are reabsorbed. Therefore 
below I will neglect the redshifts of these quanta, which could in principle take them outside of 
the resonance line. The presence of the big number of L a photons leads to an overabundance 
of the electrons (e), protons (p) and 2s, p— states of the neutral hydrogen, compared to what 
is predicted by the equilibrium Saha's formula. In turn, this delays the recombination and for 
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a given temperature the actual degree of ionization exceeds its equilibrium value given by the 
equilibrium Saha's formula. The full system of kinetic equations describing the recombination is 
rather complicated and can be solved only numerically. Here I will present the useful approximate 
treatment, which is in a very good agreement with the results of the numerical calculations. 

I neglect all highly excited hydrogen states so that the main "players left in the game" 
will be Is, 2s, 2p— states together with the electrons, protons, thermal photons, L a - and other 
energetic quanta emitted during recombination. The main processes, involving these components, 
are symbolically shown in Fig.2. They are responsible for the "converting elements" and bring 
them from one "reservoir" to the other changing the appropriate concentrations. First note that 
the recombination acts directly to the ground state do not lead to any net change in the system 
and leave the concentrations in the appropriate "reservoirs" in Fig.2 without change. 




figure 2 



Second, the thermal radiation is still very efficient and plays the dominant role in the 
ionization of the excited hydrogen atom 2s, p— states (at least, at the beginning of recombination). 
Actually to ionize the excited hydrogen atom the energy of the photon should be only one quarter 
of the binding energy Bh- The number of such photons is still bigger that the number of highly 
energetic photons emitted in recombination acts and therefore considering the ionization of the 
excited atoms one can safely ignore the distortions of the thermal radiation spectrum. On the 
contrary, these thermal quanta do not play any essential role in the transitions between Is and 
2p, s— states. The transitions Is — > 2p are mostly due to L a — quanta which at the beginning 
are present in the same number as the neutral atoms in the ground state. After the degree of 
ionization significantly drops the free electrons and the excited states are overabundant compared 
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to what one would expect according to the equilibrium Saha's formula. This is why we can neglect 
Is + 7 + 7 — > 2s transitions compared to the two-photon decay of 2s— state: 2s — > Is + 7 + 7. 
The probability of this process (W2s-»is — 8-23 sec -1 ) is much smaller than the probability of 
2p —>■ Is + L a — decay (W2 S -+i a — 4 x 10 8 sec -1 ) . Nevertheless, it plays the main role in the 
nonequilibrium recombination being, in fact, responsible for the net change of the concentrations 
of all "elements". 

The L a quanta emitted in 2p — > Is transitions are fast reabsorbed by the hydrogen atoms 
in the ground state, and these atoms go back to the u 2p— reservoir" . Therefore, the main source 
of the irreversible "leakage" from "e,p— to Is— reservoir" is the two quanta decay via 2s— levels 
and the net change of the electron concentration is mainly due to this process. All other processes 
return the "escaped" electrons very fast back to "e,p— reservoir". Hence, the rate of the overall 
decrease of the electron concentration (which is equal to the increase of the neutral atoms in the 
ground state) due to the two-photon decay of 2s— states is: 

■* = — = - 2s 2a ' (106) 

where the relative concentrations X e = n e /nt, X2 S = n2 S /nt have been introduced; here rit is the 
total number density of all neutral atoms plus electrons. I would like to stress once more that the 
equation (jl(J6|) ignores all other irreversable processes, besides of 2s — > Is + 7 + 7 decay, which 
could lead as a final outcome to the neutral hydrogen atoms in the ground state . As we will 
see later this assumption is valid until the degree of the ionization drops to rather small values. 
After that, at the end of recombination, when some other irreversable processes (in addition to 
two quanta decay) become important, I will correct the main equations to account for them. 

To express X2 S through X e let us use the quasi-equilibrium condition for "2s— reservoir". 
The rates of the reactions depicted in the Fig. 2 are very high compared to the rate of the 
expansion. Therefore the concentrations of the elements in the " intermediate reservoirs" quickly 
adjust their quasi-equilibrium values which are determined by condition that the "net flux" for 
an appropriate "reservoir" should be equal to zero. For "2s— reservoir" this condition takes the 
following form: 

(<™} ep -+ 72s U en P ~ (°)j2s->ep n 7 % 2s ~ W 2s ^i s ra 2s = 0, (107) 

where (av) are the effective rates of the appropriate reactions and ntf is the number density of 
the thermal photons. The relation between the cross-sections of the direct and inverse reactions 
can be easily found if one notes that in the state of equilibrium these reactions should compensate 
each other; hence 

where in the second equality I used the Saha's formula and took into account that the binding 
energy of 2s— state is -Bf//4. Using this relation we can express X2.S from ()107j) as 

W 2s fTm e \ 3/2 . ( V 2 
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Substituting this expression into (|1U6|) we obtain 



dt 



-w 2s 




(110) 



When the first term inside the bracket is small compared to the second one the electron and excited 
states of hydrogen atoms are in equilibrium with each other and with thermal radiation. In this 
case the last term in the equation (jl()7j) is small compared to the other terms and the relative 
concentrations of e,p and 2s— states still satisfy the appropriate Saha's relation (r.h.s. equality 
in (J108j0 . Of course, it does not mean that the ionization degree in this case is given by the 
equilibrium Saha's formula, which is derived under assumption that Is— state is also in thermal 
equilibrium with the other states. As I mentioned, the ground state drops out of equilibrium 
with the other levels soon after recombination begins and there is an overabundance of the atoms 
in the excited states compared to what one would expect according to the equilibrium Saha's 
formula 6 . 

The rate of the recombination to 2s— level is well approximated by the formula (see, for 
instance, [T2]): 



and one can easily verify that two terms inside the brackets in ()109|) becomes comparable at 
the temperature ~ 2450°fC Hence only at the temperatures higher than 2A50°K the e — p 
recombination processes are faster than the two photon decay and thermal radiation is efficient 
in keeping the chemical equilibrium between e,p and 2s— states. 

As long as the temperature drops below this value the photoionization of 2s— states be- 
comes less efficient than their two quantum decay. The thermal radiation does not play essential 
role after that and the quasi-equilibrium concentration of 2s— states is regulated by the balance 
of the recombination rate to 2s— levels and their two quanta decay rate (the second term in 
the equation (jl(J7|) can be neglected compared to the third one). In this case the second term 
inside the brackets in (jllOJ) is small compared to the first one and the rate of recombination 
due to the leakage of the electron from u e,p— reservoir" through 2s reservoir is proportional to 
(o'v) 2s n t X^ and does not depend on W 2 s ■ It is entirely determined by the rate of the re- 




combination to 2s level. At the same time 2p— states also drop out equilibrium with electrons, 
protons and thermal radiation and most of L a are destroyed in two quanta decays. As a result the 
"e,p — > 2p — > Is— channel" becomes also efficient in converting the free electrons and protons into 
the neutral hydrogen and increases the "leakage of the electrons from e,p— reservoir". Moreover, 
nearly every recombination act in one of the excited states lead to the formation of the neutral 

6 At the beginning of recombination the rate of change of the hydrogen atoms in Is— state is proportional to 
W / 2s^is^2s-This rate is much smaller than the rate of the reactions ep ^ 72s, determining the concentration of 
2s— states. 
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hydrogen atom. This effect is relevant only at the late stages of recombination and can be easily 
taken into account if we substitute in (jllUJl instead of (uv) ep ^ 2s the rate for recombination to 
all excited states, which is well approximated by the fitting formula (see, for instance, [T2"] ) 

W-^rxi«r»(*)"^ (U2) 

It is convenient to rewrite the equation (|110|) using instead of temperature and cosmological time, 
related to the temperature via (|105|) . the redshift parameter z + 1 = T/T 70 . Substituting the 
numerical values for the reaction rates and the number density n t in the obtained equation after 
elementary calculations we get: 

dX e s z \°- 3 , in4 ( 14400^-* 

^ = 15 ' 3 7Cr 7 H 14400 J + 10zex P{-—)) X e , (113) 

when we neglect unity compared to z. This equation can be easily integrated: 



XJz) =6.53x 10" 2 x 





(0.72 (z/ (1.44 x 10 4 )) - 3 + 10%exp (-1.44 x 10 4 /z)) 

One can verify that the solution X e (z) is not very sensitive to the "initial conditions" which 
could be taken at z in > z (for instance, at T ~ 3500°i^), when X e (z in ) X e (z) . The main 
contribution to the integral in this case give z < z in . At z > 900 (appropriately at the temperature 
T > 2450°i^) the first term inside the bracket in the integrand can be neglected. In this case, 
the expression ()114|) is well approximated by the f or inula [TB"]: 

X e (z) ~9.1x 10 6 ^5™ z -i exp (-1^) (H5) 
ll b h 75 \ z J 

and the overall rate of the recombination is completely determined by the rate of two quanta 
decay. It is clear from the derivation of ()114|) that this formula and, correspondingly ()115|) . are 
applicable only when the degree of ionization drops significantly below unity and the deviations 
from the full equilibrium become quite essential. As a rough criteria for the applicability of 
these formulae let us take the moment when the concentration of the neutral hydrogen reaches 
about ten percent. According to (|115j) for the realistic values of the cosmological parameters: 
a m h 2 75 ~ 0.3 and tt b h 2 7b ~ 0.03 this happens at z ~ 1220 (appropriately T ~ 3300 + 3400° K) . 
Therefore in this case the range of applicability of ()115|) is not very big, namely, 1200 > z > 900. 
However during this time when the temperature drops only from 3400°K to 2450°^ the degree 
of ionization decrease very substantially; at T ~ 2450 o fT it constitutes X e (900) ~ 2 x 10~ 2 . It is 
interesting to compare this result to the prediction of the equilibrium Saha's formula. According 
to the equilibrium Saha's formula A e (3400 o iT) ~ 10" 1 and X e (2450°f0 ~ 10" 5 , that is, at 
z ~ 900 the ionization degree exceed the equilibrium one more than in thousand times. Hence, 
the deviation from the equilibrium very essentially delay the recombination process. The other 
interesting thing is that the equilibrium ionization degree depends only on the baryon number 
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density, while in (J114)) enters also the density of the cold matter. It is not surprising, since the 
cold matter determines the rate of the cosmological expansion which is very important for kinetics 
when the deviations from equilibrium become essential. 

When the temperature drops below 2A50°K at z < 900 the approximate formula (J115|) 
is not valid anymore and we have to use (|114|) . The degree of ionization first continues to drop 
and finally freezes-out; for instance, for fl m hj 5 ~ 0.3 and Qbhj 5 — 0.03 the formula (jll4|) gives: 
X e (z = 800) ~ 5 x 10~ 3 , X e (400) ~ 7 x 10" 4 and X e (100) ~4x 10" 4 . To calculate the freeze-out 
concentration we note that the integral in (jll4|) converges for z = and is about 4 x 10 3 ; hence 



After ionization degree drops below unity the approximate results given (|114|) and (|115|) are 
in very good agreement with the numerical solutions of the kinetic equations, while the Saha's 
approximation do not reproduce the ionization behavior even roughly. 

At the beginning of recombination most of the neutral hydrogen atoms were formed as a 
result of the cascade transitions and the number of L a — photons was about the same as a number 
of hydrogen atoms. What happens with all these L a — photons afterwards? Will they survive and, 
if so, could we observe them today as an appropriately redshifted narrow line in the spectrum of 
CMB? During the recombination the number density of the L a — quanta n a is determined by the 
quasi-equilibrium condition for " L a — reservoir" 



the ionization degree squared. Thus, nearly all L a — photons which emerged at the beginning 
disappear because they are "de facto" destroyed due to the two-photons decay of 2s— states. 
Therefore there will be no sharp line in the primordial radiation spectrum. Nevertheless as a result 
of recombination this spectrum will be significantly warped in the Wien region. Unfortunately, 
the spectrum distortions lie in those part of the spectrum, where they are strongly saturated by 
the radiation from the other astrophysical sources and one can not observationally verify this 
important consequence of the hydrogen recombination. 

Finally let us find when the universe becomes transparent for the radiation. It happens 
when the typical time between the photon scattering begins to exceed the cosmological time. The 
Raleigh's cross-section for the scattering on the neutral hydrogen is negligibly small and in spite 
of their low concentration, the main role in opaqueness play the free electrons . The cross-section 
of the scattering on free electron is equal to <jt — 6.65 x 10~ 25 cm 2 and the equation defining the 
moment when the radiation completely decouples form the matter takes the form: 
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W 2 p^i s n 2 p = (cr a )n a n ls . 
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Since n ls « n t and n 2p oc X 2 



we see that the number of these quanta drops proportionally to 



1 



(118) 



a T n t X e 



35 



Xf c ~ 40-^5 ( ZkY (i .| «.)) 



This equation can be rewritten as 

3/2 

^6^75 \TdecJ 

By "try-out" one can easily check that the decoupling happens at T^ ec ~ 2500°/^ ( the corre- 
sponding redshift Zd ec ~ 900) irrespective how big are the values of the cosmological parameters. 
If Vt m h^ b ~ 0.3 and r2 fe /i| 5 ~ 0.03 the ionization degree at this moment is about 7 2 x 10~ 2 . 



B Asymptotic behavior of the transfer functions 

The resulting fluctuations of the background radiation depend on the gravitational potential $ 
and the radiation energy density fluctuations <5 7 = fe 7 /e 7 at the moment of recombination. To de- 
termine these quantities we have to study the gravitational instability in two component medium 
consisting of the coupled baryon-radiation plasma and the cold dark matter. Because these com- 
ponents interact only gravitationally their energy-momentum tensors conserve separately. In the 
cosmological conditions the shear viscosity can not be neglected for the baryon-radiation plasma 
and leads to the dissipation of perturbations in small scales (Silk damping). For imperfect fluid 
with the energy density e and the pressure p one can use the energy-momentum tensor given in 8 
Then one can easily find that in a homogeneous universe with small perturbations described 
by the metric (0) the conservation laws Tg. a = in the first order in perturbations reduce to 

Se' + m (5e + 5p) - 3 (e + p) $' + a (e + p) u { >t = 0. (120) 

\ (a 5 (e + p) Vb - -riAu* >t + ASp + (e + p) A$ = 0. (121) 

where 8e, Sp are, respectively, the perturbations of the energy density and pressure ; u % is the 
peculiar 3- velocity and i) is the shear viscosity coefficient. Note that the first equations which 
follows from Tq. q = does not contain the shear viscosity. The second equation was obtained by 
taking the divergence of the equations T? a = 0. As it was already noted, these two equations are 
separately valid for the dark matter and for the baryon-radiation plasma. 

Dark matter. For dark matter, the pressure p and the shear viscosity 77 are both equal to 
zero. Taking into account that e^a 3 = const we obtain from ()120j) that the fractional perturbations 
in the energy of dark matter component 5d = 5ed/£d satisfy the equation 

(5 d -3$)' + au\ i = 0. (122) 

7 It is rather interesting to note that this time coincides with the moment when e,p and 2s— levels come out of 
equilibrium and the approximate formula (|115J) becomes inapplicable. 

8 I will neglect the heat conduction since it does not change substantially the Silk damping scale. 
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If we express u\ in terms of 5^ and $ and substitute it into (|121|) then the resulting equation 
takes the following form: 

(a(5 d -3$)')'-aA$ = 0. (123) 

Radiation-baryon plasma. The baryons and radiation are tightly coupled before recom- 
bination and, therefore, generically only the sum of their energy-momentum tensors satisfies the 
conservation laws (J12U)) and (|121|) . Nevertheless, in particular case when the baryons are nonrel- 
ativistic, the equation (|12(J|) is still valid separately by baryons and radiation, because the energy 
conservation law for the baryons, T^ a = 0, reduces in this case to the conservation law for the 
total baryon number. (Of course this is not true for ()121Jl ) since baryons and radiation "move 
together" and there exist momentum exchange between these components.) Hence, the fractional 
density fluctuations in baryons, 5b = 5Eb/£b, satisfy the equation similar to (|122|) : 

{S b - 3$)' + au\i = 0. (124) 

As it follows from (J12U)) the corresponding equation for the perturbations in the radiation com- 
ponent, 5 7 = Se-Je^ takes the form 

4 

(<5 7 -4$)' + -au\i = 0. (125) 

Since the photons and baryons are tightly coupled their velocities are the same. Therefore mul- 
tiplying ()125J) by 3/4 and subtracting equation ()124|) we obtain 

5s 3 

— = -5 7 — 5b = const (126) 

where 5s /s are the fractional entropy fluctuations in the baryon-radiation plasma. For adiabatic 
perturbations, 5s = and, therefore, we have 



<V, = !* r 1127) 



If we express u\ in terms of <5 7 and $ from (|125|) and substitute into (|121|) . we obtain 

'4)'_^ H _ A , 7= jL A$+ (if:)'_^, (128) 

c~ / e 7 a ' 3c; \ c% J e 7 a 

where A is the Laplacian and c 2 s is the squared speed of sound in the baryon-radiation plasma, 
which is equal to: 

sj, = sp^ = U 1 + 3s 1 y\ 

5e 5e 1 + 5E b 3 \ Ae~J v ' 

Without taking into account the polarization effects the shear viscosity coefficient entering ()128j) 
is given by[7j: 

4 

f] = -^ e -i T -i ( 13 °) 
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where r 7 is the mean free time for the photons. 

Thus we derived two perturbation equations ()123|) and (|128jh which being supplemented 
by 0-0 component of the Einstein equations [3] 

A$ - m& - 3H 2 $ = AnGa 2 (e d 6 d + ^ e -r S ij ( 131 ) 

form a closed system of equations for three unknown variables 8^, 5 1 and $ (we used (I127|) to 
express 6b in terms of <5 7 ). 

^From f)125j) it follows the useful relation for only the radiation contribution to the diver- 
gence of — i components of the energy-momentum tensor, 

Tl i = ^e 1 u Q u i ti ={^-5 1 )'e 1 (132) 

which is used in (J2UJ) 

Longwave perturbations: (k <C Vr 1 ) The behavior of perturbations strongly depends on 
how big is their scales compared to the horizon. First I consider the long wavelength perturbations 
with ki] r <C 1 (k is comoving wavenumber) which cross the horizon only after recombination. 
Knowing the gravitational potential we can easily find £ 7 . In fact for this longwave perturbations 
one can neglect the velocity term in the equation ()125j) . which after that can be easily integrated 
with the result 

5 7 - 4$ = C, (133) 

where C is the constant of integration. To determine C, we note that, during the radiation 
dominated epoch, the gravitational potential is mostly due to the fluctuations in the radiation 
component and does not change on supercurvature scales. At early times, <5 7 ~ —2$ (77 <C r) eq ) = 
—2$° (see [3]) ; hence C = —6$°. After equality, when the dark matter overtakes the radiation, 
the gravitational potential $ changes its value by factor of 9/10 and then remains constant, that 
is, $ (77 ^> r] eq ) = (9/10) $°. Therefore, if cold dark matter dominates at recombination, it follows 
from (fH31) that 

5 1 ( Vr ) = -6$° + 4$ ( Vr ) = -h (r) r ) . (134) 

One arrives at the same conclusion by noting that, for the adiabatic perturbations, <5 7 = 45^/3 
and 5d ~ —2$ (r] r ) at recombination. 

Intermediate scales (r/" 1 < k < r/" 1 ) Next I consider the scales which enter horizon in 
between the equality (i] eq ) and recombination (r] r ). The perturbations which enter horizon within 
this rather short time interval are especially interesting since they are responsible for the first few 
acoustic peaks in the CMB spectrum. Unfortunately, for the realistic values of the cosmological 
parameters the solution for these perturbations cannot be found analytically with needed accuracy 
because in the realistic models the condition 7] eq rj r is not satisfied. Nevertheless to gain an 
intuition about the behavior of perturbations, it is very useful to consider the models where r\ eq 
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r\ r and derive the appropriate asymptotic expressions for the perturbations with n~ x <C k <C n~ x . 
To simplify the consideration I also assume that the contribution of baryons to the gravitational 
potential is negligible compared to the contribution of the cold dark matter. 

In general, there exist four instability modes in two component medium. The set of 
equations is rather complicated and they can not be solved analytically without making further 
assumptions. However in our case, the problem can be simplified if we note that if the perturbation 
enters horizon sufficiently late after equality (77 rj eq ), the appropriate gravitational potential, 
which is mainly due to the perturbations in the cold dark matter component, remains unchanged 
and stays constant afterwards (if) = const) [3]. The baryons do not contribute much to the 
gravitational potential, however they can still significantly influence the speed of sound after 
equality. 

Under assumption we have made, the gravitational potential $ can be considered as an 
external source in equation (|128j) . Therefore, the general solution of this equation is given by 
the sum of a general solution of homogeneous equation (with $ = 0) and a particular solution of 
(I128p . Introducing the variable x, defined by dx = c 2 dn, and taking into account that the time 
derivatives of the potential ()128|) are equal to zero ($ = const), we reduce the equation (|128J) to 

d S'y 4t*^y d/Sty 1 ^^(^ ^ /\ ^135^ 
dx 2 5a dx c 2 s 7 3cf 

where the second term is due to the viscosity. If the speed of sound is slowly varying, this equation 
has an obvious approximate solution 



5, ~ (136) 



The general solution of the homogeneous equation ()135|) can be obtained in the WKB approxima- 
tion. Let us consider the plane wave perturbation with the comoving wavenumber k. Introducing 



instead of 8 1 the new variable 



y = 5 7 exp ( -k 2 I —dx ) , (137) 




5 

we find from ()135j) that it satisfies the equation 

sH^P?) i»= a « i38) 

For the perturbations with the scale (X p h ~ a/k) bigger than the mean free path of the photons 9 
(~ r 7 ), the second term inside the brackets is negligible. The third term which is about r^/arj ~ 
T 7 /t 1 can be also skipped. Therefore the WKB solution for y is 

ClCO Jk [ —)+C 2 sin(k [ —)). (139) 



c, / \ . c 



In fact, the imperfect fluid approximation can be used only in this case 
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Returning back to <5 7 (see (|137)0 and combining this solution with (|136|) . we obtain 

5 y ~ + V^s (c x cos (k J c s d V ^j + C 2 sin (k J c s d V j ) e^^ 2 . (140) 

Here we have introduced the dissipation scale characterized by the comoving wavenumber: 

k D (rj)^(^j\^dr^j V2 . (141) 

In the limit of constant speed of sound and vanishing viscosity the solution ()140|) is exact and 
valid also in the limit k — > 0. 

^From ()140j) . it is clear that the viscosity efficiently damps the perturbations on comoving 
scales A < 1/ko- Using the formula (|141j) with c\ = 1/3 and, assuming instantaneous recombi- 
nation we obtain the following estimate for the dissipation scale 

(kovr)- 1 ~ 0.6 (O^ 2 ) 174 {n b h 2 y 1/2 z; 3/A - (142) 



The constants of integration C\ and C 2 in ()140|) can be determined if we note that 
at earlier stages when the speed of sound does not change too much the solution ()140|) is also 
valid when the scale of perturbation still exceeds the horizon scale. As we have found before the 
amplitude of the longwave perturbations ikn <C 1) is equal to <5 7 ~ — 8$fc/3 = const at n 3> n eg . 
Assuming that at the moment when the perturbation enters the horizon the speed of sound is 
still not very different from 1/a/3 we find that C\ = 4/3 3 / 4 and C 2 = 0; hence 

^l) (143) 







for k <C r}~ . Here we took into account that = 9$yl0 and expressed the result in terms 
of the initial gravitational potential on superhorizon scales before equality The result ()143j) 
coincides with the result obtained by S. Weinberg [2| in synchronous coordinate system. 

Shortwave perturbations (k ^> 77" ) Finally I consider the perturbations which enter 
the curvature scale before equality. At n <C r] eq , the radiation dominates and in this case the 
appropriate expressions for $ and 5 7 were derived, for instance, in Neglecting the decaying 
mode we find that after perturbation entered the horizon, that is, at k~ x <C i] <C n eq 

5 1 ~ 6$£ cos ( krj/Vs) , $ fc (n) ~ — ^% cos ( kn/Vt] . (144) 
v ' {kn) v ' 

The dissipation, which becomes important only before recombination, can be treated similar to 
how it was done above. Therefore I neglect the dissipation term here and restore the damping 
factors only in the final expressions. 
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After inhomogeneity entered the horizon, the cold dark matter starts "to slide" with 
respect to the radiation. To get an idea about the behavior of the inhomogeneities in the cold 
dark matter component itself we can use the equation ([123)1 . which after integration becomes 

5 d = 3$ + J ^- J aA^drj". (145) 

Note that this is an exact relation which is always valid for any k. During the radiation-dominated 
epoch, the main contribution to the gravitational potential is due to the radiation, and, therefore, 
the gravitational potential in the equation ()145|) can be treated as an external source. We can 
fix the constant of integration in (|145|) substituting there the exact solution for the radiation 
dominated universe (see (5.45)-(5.46) in [3|) and noting that, at earlier times on superhorizon 
scales one has to match the well-known result for the longwave perturbations: 5d — 3<5 7 /4 ~ 
— 3$°/2. As a result we obtain that after entering the horizon, but before equality 



5 d ~ -9 (c-~ + In (krj/Vt) + O ((fc^ 1 )) $° fc , (146) 

where C = 0.577... is the Euler constant. That is the perturbations in the cold matter component 
are "frozen" (they grow only logarithmically) 

^From (jl3ip it is easy to see that before equality the contribution of dark matter perturba- 
tions to the gravitational potential is suppressed by the factor £d/e 7 compared to the contribution 
from the radiation component. At equality, the dark matter begins to dominate and the density 
perturbation 5d grow as oc r] 2 (see, |3j). As a result, the appropriate gravitational potential "freeze 
out" at the value 



, , . AnGa 2 € . 

®k [V > Veq) p d 



O( l)^|)$o (147) 



and stays constant until the recombination. One can get the exact numerical coefficients in 
this formula in the following way. For shortwave perturbations, the time derivatives of the 
gravitational potential in (|123|) . ()131|) can be neglected compared to the spatial derivatives. Then 
from these equations it follows that 

{aS' d y - AnGa 3 (s d 8 d + ^j^t) = °- ( 148 ) 



The second term here induces the corrections to the solution ()146|) which become significant only 
near equality. These corrections are mostly due to e d 8d— term. The term oc e 7 5 7 does not lead 
to essential corrections to the solution (|146|) before equality and it is also negligible compared to 
£ d^d~ term after equality; hence it can be skipped in ()148|) . As a result the obtained equation can 
be rewritten in the following form 
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where x = a/a eq . The general solution of this equation is (see, for instance, jjj) 



At x <C 1 it should coincide with 



1 + r 



VI + x 
. Comparing ()150|) with 



In • ' - — 



Ci ~ -9 In 



7T 



c 



$2 



c 2 



at i C 1 we find 



9<&° 



(150) 



(151) 



where 77* = rj eg / (V2 — l) . During the matter dominated epoch (x 1), the second term in ()150j) 
corresponds to the decaying mode. Neglecting this mode, assuming that the baryon contribution 
to the potential is negligible compared to the dark matter and using the relation between the 
gravitational potential and 8d (see (I131j) ) one finally gets 



$fc (?7 > Veg) 



In (0.15% e? ) $0 



(0.27% 



(152) 



eq 1 



in agreement with The fluctuations in the radiation b~ 1 after equality continue to behave as 
sound waves in the external gravitational potential given by ()152j) . Therefore, they are described 
by ()140J) . where we have to substitute the potential ()152j) instead of The constant of inte- 
gration can be fixed by comparing the oscillating part of this solution to the result in ()144j) at 
r] ~ r] eq . Then, we find that at r\ r\ eq 



5 7 ~ 



4 In (0.15%, 



eq I 



3c2 (0.27%, 



eq I 



+ 3 5/4 \/4c^cos i^k J c s dT]^j e 



-(k/k D f 



(153) 



for k 3> rj eq . We have restored here the Silk damping factor. During the radiation dominated 



epoch, the damping scale, which is proportional to the photon mean free path, is very small. 
However, it increases just before the recombination and therefore the oscillating contribution to 
5 7 is exponentially suppressed on small scales. 
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